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The sensitivity of all-sky searches for gravitational-wave pulsars is primarily limited by the finite 
availability of computing resources. Semicoherent searches are a widely-used method of maximizing 
sensitivity to gravitational-wave pulsars at fixed computing cost: the data from a gravitational-wave 
detector are partitioned into a number of segments, each segment is coherently analyzed, and the 
analysis results from each segment are summed together. The generation of template banks for 
the coherent analysis of each segment, and for the summation, requires knowledge of the metrics 
associated with the coherent and semicoherent parameter spaces respectively. We present a useful 
approximation to the semicoherent parameter-space metric, analogous to that presented in Wette 
and Prix [Phys. Rev. D 88, 123005 (2013)] for the coherent metric. The new semicoherent metric is 
compared to previous work in Pletsch [Phys. Rev. D 82 , 042002 (2010)], and Brady and Creighton 
[Phys. Rev. D 61, 082001 (2000)]. We find that semicoherent all-sky searches require orders of 
magnitude more templates than previously predicted. 

PACS numbers: 04.80.Nn, 95.55.Ym, 95.75.Pq, 97.60.Jd 


I. INTRODUCTION 


for post-processing promising detection candidates [e.g. 


Gravitational-wave pulsars [I] are rotating neutron 
stars which are hypothesized, though not yet observed, 
to emit gravitational radiation in the form of long-lived, 
continuous, quasisinusoidal signals. The maximum am¬ 
plitude of a gravitational-wave pulsar signal has been 
most recently studied in [5], and the evolution of the 
population of Galactic neutron stars emitting at an as¬ 
sumed amplitude has been modeled in [311]. The sen¬ 
sitivity of the interferometric gravitational-wave obser¬ 
vatories LIGO 13 m and Virgo are expected to 

improve significantly within this decade, while a third 
large-scale observatory, KAGRA [9|, is under construc¬ 
tion. Nevertheless, it remains likely that gravitational- 
wave pulsars will remain at best marginally detectable 
for the foreseeable future. 

The detection of a gravitational-wave pulsar, in the ab¬ 
sence of an observed electromagnetic counterpart [cf. [13 , 
is a significant data-analysis challenge, requiring both 
large-scale computing resources and highly-optimized 
data analysis algorithms. The former is provided by 
the distributed computing project Einstein@Home El, 
which is dedicated to searching for both gravitational- 
wave and electromagnetic pulsars; to date, the project 
has discovered many new radio and gamma-ray pul¬ 
sars [e.g. m II3. Research into the latter has led to 
advances in the understanding of optimal and/or ro¬ 
bust detection statistics [e.g. ITdHT^ . the construction 
of matched-filtering template banks [e.g.|IlM3, hierar¬ 
chical search methods [e.g. ES1 - E5] and their optimiza¬ 
tion given limited computing resources [e.g. 123 EO], the 
estimation of search sensitivity [e.g. ED, and pipelines 
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This paper builds on recent work on all-sky, broad¬ 
band, coherent searches for isolated gravitational-wave 
pulsars. (Similar studies of the parameter-space metric 
for gravitational-wave pulsars in binary systems can be 
found in e.g. [311|3S|-) In (23 (hereafter Paper I), a use¬ 
ful approximation to the parameter-space metric is devel¬ 
oped, which is important for the correct construction of 
banks of matched-filtering templates. Previous work [e.g. 
123113 had assumed restrictions on the time-span of data 
which could be coherently analyzed; these restrictions 
were relaxed in Paper I In [22] (hereafter Paper III, a 
practical implementation of an optimal template bank 
using lattice template placement is presented. The num¬ 
ber of templates in the bank was found to be in good 
agreement with the prediction of m- 


This paper presents a useful approximation to the 
parameter-space metric for all-sky semicoherent searches 
for gravitational-wave pulsars, building on the work 
of IPaper I[ In a semicoherent search, the data are divided 
into a number of segments, each of which are coherently 
searched using a per-segment coherent template bank] the 
search results from each segment are then summed to¬ 
gether using a semicoherent template bank [38] • Each 
semicoherent template preserves a consistent (instanta¬ 
neous) signal frequency fit) between each segment, but 
not necessarily a consistent signal amplitude. While 
semicoherent searches are less sensitive than a single¬ 
segment fully-coherent search of the same data, they are 
also computationally cheaper, and are therefore able to 
analyze data spanning years or more, considerably longer 
than would be possible using coherent methods alone. 

Section [^ of this paper reviews pertinent background 
information. Section IIIII derives the new semicoherent 
parameter-space metric, and validates it using numeri¬ 
cal simulations. Section jlVj compares the new metric to 
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Quantity 

Symbol 

Number of segments 

N 

Segment index 

1 

Phase evolution parameters of signal 

A" 

Difference with respect to signal parameters 


Time-span of coherently-analyzed data 

f 

Parameters of coherent template 

Xt 

Parameters of nearest coherent template 

Xe{---) 

Goherent J^-statistic 

2Tt 

Goherent signal-to-noise ratio 

SNRe 

Goherent noncentrality parameter 

P? 

Goherent 7^-statistic metric 

g° 

Goherent J^-statistic mismatch 


Goherent phase metric 

~d) 

SI 

Goherent reduced supersky metric 

Se 

Goherent reduced supersky metric mismatch 

Jli 


Time-span of semicoherently-analyzed data 
Parameters of semicoherent template 
Difference with respect to semicoherent template 
Semicoherent 7^-statistic 
Semicoherent signal-to-noise ratio 
Semicoherent noncentrality parameter 
Semicoherent 7^-statistic metric 
Semicoherent J^-statistic mismatch 
Semicoherent phase metric 
Semicoherent reduced supersky metric 
Reduced supersky mismatch without interpolation 
Reduced supersky mismatch with interpolation 
Semicoherent global correlation metric 
Global correlation mismatch without interpolation 
Global correlation mismatch with interpolation 


T 

A 

A 

2T 

snr"^ 

^2 


!</> 

§ 

fJ-ni 


TABLE I. Summary of notation used in this paper. The 
decoration of signal, coherent, and semicoherent quantities 
with andT respectively follows the convention of 130 ]. 


previous work [211 [21] , and Section |v] summarizes the 
conclusions of the paper. 


II. BACKGROUND 


This section briefly reviews basics of gravitational- 
wave pulsar searches (Section IIA I, a nd intro duces 
the coherent parameter-space metric of [Paper I (Sec¬ 


tion IIB). It then gives a summary of the theory and 
notation of semicoherent searches (Section II C[ ); for a 
similar review see |30j . The notation introduced in this 
and subsequent sections is summarized in Table [^ 


A. Coherent gravitational-wave pulsar searches 


The signal from a gravitational-wave pulsar observed 
in a detector at time t is modeled m by the function 
h{t, A., A). The four parameters A control the amplitude 
modulation of the signal, and the vector of parameters 
A control its phase evolution. For all-sky searches for 
isolated pulsars, A includes the sky position of the pulsar, 
its initial frequency / = = /(to) at some reference 

time to, and frequency derivatives (or spindowns) = 


A coherent search using the J--statistic [B M 
matched-filters data from a gravitational-wave detector 
against h(t, A, X), and analytically maximizes the log- 
likelihood j;atio over the A. The resultant detection 
statistic 2J'((X) is thus a function only of A. (The sub¬ 
script -i here is the segment index, which will become 
relevant for semicoherent searches in Section II C| ) To 


further maximize 2J-i(X) over the parameter space V of 
possible X GV, a finite set of parameters {A^} C 7^ - the 
coherent template bank - is generated, and 2iFi{Xg) is 
computed for all Xi G {A^}. The most promising pulsar 
candidate, assuming Gaussian detector noise [cf. [18], will 
be the maximum of 2J^^(A^) over {A^}. 

We adopt the following definition for the signal-to- 
noise ratio: 


2 _ E[2J^ near signal] — E[2J^ in noise] 

SNR =- - - 

^jYav\^F in noise]/8 


( 1 ) 


where E[-] denotes the expectation and Var[-] the vari¬ 
ance. (The rationale for the normalization of the vari¬ 
ance by 8 can be seen in Eqs. and ( [T4| below.) In 
noise, 2iFe{Xi) is a central statistic with 4 degrees 
of freedom; near a signal, it is a noncentral statistic 
with 4 degrees of freedom and noncentrality parameter 
'^({A, A"; A^), where A^ denotes the phase parameters of 
the matched-hlter template and A® the true phase pa¬ 
rameters of the signal. It follows that the signal-to-noise 
ratio of the coherent J^-statistic is 


SNrI{A, X^;Xe) 


[4 -I- pj{A, A**; Af)] 

vW8 

p?(AA^A,). 


[4] 


( 2 ) 


In the gravitational-wave pulsar literature, pf is often 
also referred to as the signal-to-noise ratio; however, as 
noted in |3S] , this is true only for a fully-coherent search, 
as will be seen in Section fll Cl 

Since the true phase parameters of a pulsar, denoted 
A®, will never exactly match one of the {A^}, some loss of 
signal-to-noise ratio is inevitable. The mismatch quan¬ 
tifies the loss in signal-to-noise ratio SNR^(^, A**; A^) at 
a template A^ ^ X^, relative to that of a perfect match 
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SNR^ (AA^A^) ; it is defined to be 


= (3) 

SNR^(A A"; A") 

Substitution of Eq. ([^ gives the mismatch alternatively 
in terms of the noncentrality parameters: 


/20(AA^A,) = 1- 


P?(AA^A,) 

p?(Aa^a^) 


(4) 


The definition of mismatch in terms of loss of signal-to- 
noise ratio is consistent with the data analysis of searches 
for signals from compact binary coalescence [e.g. 
definition in terms of loss of noncentrality parameter is 
more common in gravitational-wave pulsar data analy¬ 
sis [e.g. [30]. In this paper, by taking Eq. (|^ as the pri¬ 
mary definition of mismatch, we aim to highlight the dif¬ 
ference between signal-to-noise ratio and noncentrality 
parameter for semicoherent searches (see Section IIC). 
Ultimately, it will be seen that the two definitions of mis¬ 
match given by Eqs. (§ and @ are equivalent. 

Efficient template bank construction seeks to minimize 
the average Jij, over all possible A**, at fixed number of 
templates A/f = |{A^}|; this can be ac hieved using lat¬ 
tice template placement (see Paper II |. To aid in con¬ 


structing the template bank, the parameter-space metric 
is introduced Enmoj. A second-order Taylor expansion 
of Eq. Q with respect to small parameter differences 
A^Xi = Af — A® gives 


^^(A, A^ A^) « A^A, • i°(A A") • A^A, 


(5) 


where • denotes the vector/matrix inner product. The 
parameter-space metric is represented by the matrix of 
components 


i°(AA^) 


-1 ap^(AA-;A) 

2p2(XA';A'^) d\ 


( 6 ) 


and is a function of A. and A^. 

Lattice template placement, however, requires Eq. (§ 
to be constant with respect to ^ and A useful ap¬ 
proximation which partially satisfies this requirement is 


M?(A A^ a,) « A^A, • g,^A^) • A^A,; (7) 

the phase metric gf(A®) approximates g£(^, A®), and is 
given by 


gf{X) = (5a^ 0 - {dx(l))t ® (9a(/>)Ja=a= , (8) 

where dx is the vector of derivatives with respect to A, ® 
denotes the vector outer product, and the {■)t operator 
time-averages over the segment time-span T. The pulsar 
phase (/)(t. A) is given by (e.g. [Paper Ij) 






s=0 


jt - to)-+i 

(s + 1)! 


r{t) 


./max ; 


(9) 


where r{t) is the detector position relative to the So¬ 
lar System barycenter, n is the pulsar sky position, and 
/max is ^ constant chosen conservatively to be the max¬ 
imum of f{t) = X)/dt over T. It follows from 

Eq. ([^ that gf(A®) will also be constant with respect 
to any phase parameter A® in which A) is linear. 
As seen from Eq. (|^, such parameters include the fre¬ 
quency / and spindowns but not the sky position 
n = (cosacosi5, sinacosd, sinJ) when parameterized in 
terms of right ascension a and declination b. 


B. Coherent supersky parameter-space metric 


Paper I presented a close approximation to g|’(A®) 
rich is constant with respect to all phase parameters 
A, including sky position: the reduced supersky metric g^ 
(denoted grss in|Paper I). It is derived from gf(X^) by 
the following procedure, simplified from [Paper I| 

(i) gf(X^) is computed using the three components of 
n in equatorial coordinates as sky coordinates; this 
metric is the unconstrained supersky metric g^ (de¬ 
noted gss in Paper I). Since (/(t. A) is linear in n. 


this yields a metric which is constant with respect to 
all phase parameters; however, the two-dimensional 
sky is now embedded in three dimensions, which is 
undesirable for template placement. 


(ii) A linear coordinate transformation is performed on 
g^ which exploits the near-degeneracy between the 
signal phase’s dependence on the Earth’s orbital 
motion, and its dependence on the pulsar’s fre¬ 
quency and spindowns. The transformation also 
removes from g“ any correlation between sky and 
frequency/spindown coordinates. 

(iii) The eigenvectors and eigenvalues of the 3x3 block of 
the transformed g” pertaining to the sky parameters 
n are computed. The eigenvector corresponding to 
the smallest eigenvalue is identified; by definition, 
offsets A®Af parallel to this eigenvector will change 
the mismatch A^Afg^-A^A^ by the smallest amount 
possible. 

(iv) The transformed metric g“ is projected onto the 

hyperplane in the coordinate space of n perpen¬ 
dicular to the eigenvector identified in step ( pli| . 
This reduces the number of sky parameters back 
to two, and yields the metric gg. The phase pa¬ 
rameters associated with gg are the sky parame¬ 
ters {T^i,ng,g), and the frequency/spindowns Vg^'^ = 
jis) _|_ where the A^g are vectors computed 

during step (§■ 

We denote the linear transformation of the components 
of gg described in step © by the operator T, and the 
projection described in step ( jf^ by the operator ^7, i.e. 

gg = jT[ri] ■ ( 10 ) 


27r 
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FIG. 1. Relationship between the signal parameters A®, its 
nearest semicoherent template A, and its N nearest coherent 
templates Ai, A 2 ,..., A^,..., following Fig. 1 of (30| . 


[Paper I] demonstrated, using numerical simulations, 
that is a good approximation to g^(A®) for segment 
time-spans T > 1 day, maximum mismatches Jif ^ 0.6, 
and spindowns up to 


C. Semicoherent gravitational-wave pulsar searches 

This section describes a common semicoherent search 
method for gravitational-wave pulsars, referred to as “in¬ 
terpolating StackSlide” in |5D] . 

The gravitational-wave detector data are partitioned 
in time into N segments, labeled with the index £ = 
1, 2,..., For simplicity we assume that the time-spans 
of all segments are identical and equal to T ; the mid-time 
of each segment is denoted ti. In each segment £, the 
coherent J^-statistic 2J^^(A^) is calculated on a coherent 
template bank {A^} G V which is constructed using the 
coherent metric at a common reference time to- As 
the metrics g^ are generally different for each segment, 
so too will be the template banks {A^}. 

The semicoherent A-statistic is computed on a distinct, 
semicoherent template bank {A} G V. For each semico¬ 
herent template A, the coherent template in each segment 
which is nearest to A, as measured by the respective co¬ 
herent metric g^, is determined; these N templates are 
denoted The metric mismatch between A and 

A£(A) is denoted 

A,(A)) = AA,(A) • g, • AA,(A), (11) 

where AAf(A) = Af(A) — A; see Fig. Finally, the 
N coherent J^-statistics computed at the A^(A), denoted 


2J-,(A,(A)), are retrieved and summed to give 

N 

2.F(A,(A)) = ;^2.F,(Af(A)) . (12) 

i=i 

In noise, 2J^(A£(A)) is a central statistic with iN 
degrees of freedom; near a signal, it is a noncentral x^ 
statistic with 4N degrees of freedom and noncentrality 
parameter 


N 


p2(A;A,(A))=^^(A;A,(A)), (13) 


where ^(A®;A^(A)) are the noncentrality parameters of 
the N coherent J^-statistics. (Henceforth we suppress the 
dependencies on the A..) By Eq. Q, the signal-to-noise 
ratio of the semicoherent J^-statistic is therefore 


SNR^(A'^;A^(A)) 


[4A + p^(A;A^(A))] 

vWP 


[4A] 


(14) 


Note that Eq. (141 is a factor of vV less than the signal- 


to-noise ratio of a fully-coherent search using the same 
amount of data NT. 

The definition of the semicoherent mismatch follows 
from Eq. 


fJ- 


(V;A,(A)) = 1- 


SNR (A^A^(A)) 
S1?R^(A";A") 


(15) 


where SNR (A®; A**) is the signal-to-noise ratio at per¬ 
fect match in both coherent and semicoherent templates. 
Substituting Eq. ( |T4| demonstrates that the semicoher¬ 
ent mismatch can equivalently be defined as the loss in 
noncentrality parameter, c.f. Eq. (18) of [30] : 


/20(A;A,(A)) =1- 


p^(V;A,(A)) 

^(V;V) 


(16) 


where (A**; A®) is the noncentrality parameter at perfect 
match in both coherent and semicoherent templates. It 
is assumed that the coherent noncentrality parameter at 
perfect match is constant over all segments, which implies 


f{N-,\^)=Npj{N;N). (17) 

Substituting Eqs. ( |l^ , and 0 into Eq. yields 


/I°(A;A,(A)) =1- 


1 pf (A ; A^(A)) 
iV^ P?(V;V) 


1 

N 


N 

^/I0(V;A,(A)) 


e=i 


(18) 

(19) 
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where (A**; A^(A)) is the coherent mismatch in the ^th 

segment between the nearest template A^(A) and the sig¬ 
nal A®. Further substitution of Eq. 0 results in 


N 


^°(A^A,(A)) « -5]A^A,(A).g,^A^).A^A,(A), (20) 


where (see Fig. 


A'^A^(A) = A^(A) - A*^ 

= (A,(A) - A) + (A - A^) 
= AA^(A) A"A. 


( 21 ) 

( 22 ) 

(23) 


Finally, substituting Eq. (23) into Eq. (20) gives 

N 




(A^A,(A))«l^A^A.i,^A^).A^A 


N 


1^1 

N 


-^AA,(A).g,^A^).AA,(A) (24) 


N 

2 

N 


e=i 

N 


^AA,(A)-g,^A^)-A^A. 


e=i 


The first term on the right-hand side of Eq. (24) is 


the mismatch between the signal A and the semicoher- 
ent template A in the ideal case of no interpolation, i.e. 
supposing AA^(A) = 0 for all It may be written as 




N 


N 

1 

N 


^=1 

N 


^A«A.g,^A^).A^A 


^=1 


A"A-g'^(A") • A"A, 


where 


N 




(25) 

(26) 

(27) 

(28) 


e=i 


is the semicoherent parameter-space metric; cf. Eq. (2.9) 
of [221. The second term on the right-hand side of 


Eq. (24) is the mismatch purely from interpolation, 
1 ^ 


)x°i(A;A,(A))^-^)I°(A;A,(A)) 


(29) 


^=1 

N 


-5]AA,(A).g,^A^).AA,(A). (30) 


1=1 


The average of the third term on the right-hand side 
of Eq. (24) will tend to zero when N is large, and when 


/iT is averaged over many different signals A [30] . This 


is because the term is linear in the offsets AA£(A) and 
A®A, which are presumed to be equally distributed about 
zero. Likewise, under the same conditions the average of 
Eq. (29) will tend to (^) , the average coherent mismatch. 
The average of Eq. (15) will therefore tend to 


(/2«)^(/2°i) + (p(°). 


(31) 


It follows from Eqs. ( |31| ) and (25) that the average semi¬ 
coherent mismatch IJP) is minimized if A S {A} is chosen 
to be nearest to the signal A®, as measured by the semi¬ 
coherent metric , thus minimizing (/x^i)- 

Therefore, the semicoherent template bank {A} should 
be constructed similarly to the coherent template banks 
{Af}, using g'^ to minimize /x®;, and hence (/x°), at fixed 
number of semicoherent templates M = |{A}|. This may 
be ach ieved using lattice template placement (see Pa- 
, but only if g"^ is constant with respect to A. and 
le construction of such a metric, starting from the 


per II 

semicoherent phase metric g'^(A®) = is 

the subject of the next section. 


III. SEMICOHERENT SUPERSKY 
PARAMETER-SPACE METRIC 

This section introduces a close approximation to 
g‘^(A'*) which is constant with respect to all phase pa¬ 
rameters A: the semicoherent reduced supersky metric 
g. Its derivation, presented in Section III A closely fol¬ 
lows that of the coherent metric g^, which is described 
in Section |IIB| and [Paper I| The new metric is validated 
using numerical simulations presented in Section [III B| 


A. Derivation 

First, gf (A®) are computed for all segments, at a com¬ 
mon reference time to, using the components of n as sky 
coordinates; these metrics are the coherent unconstrained 
supersky metrics g^. Second, the g^ are averaged to give 
g'^ = semicoherent unconstrained su¬ 

persky metric. Finally, the procedure described in Sec¬ 
tion IIB and Paper I[ denoted by the operators JT, is 
applied to g'^ to give the semicoherent supersky metric 


g = 


JT 

N 


■ N 

.e=i 


(32) 


The phase parameters associated with g are the sky pa¬ 
rameters (na,ni,), and the frequency/spindowns = 
y(s) _|_ where the are vectors computed when 

applying T. 

Figure|^illustrates this process using an example setup 
of A = 5 segments, each spanning T = 4 days, with seg¬ 
ment mid-times tg distributed evenly over a total time- 
span T = 364 days. The effect of the changing segment 
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FIG. 2. Metric ellipsoids of the coherent and semicoherent supersky metrics, at 5 example sky points, (a) ■ (e) Metric 
ellipsoids (red) of 5 coherent unconstrained supersky metrics g“, in equatorial coordinates {x,y,z), for /Imax = 30, T = 4 days, 
/max = 1 kHz, to = UTC 2015-01-01 00:00:00, and ti — to = -180,-90,0,-1-90,-1-180 days respectively. Their intersection 
with the sky sphere |n| = 1 (blue) reproduces the coherent physical sky metric (yellow). |(f)| Metric ellipsoids (red) of the 
semicoherent unconstrained supersky metric g“ = Sr/3, for Pmax = 1500. Its intersection with the sky sphere (blue 

reproduces the semicoherent physical sky metric (yellow). 


(g) 


coordinates (a, b, c) given by its sky eigenvectors, for /2 
semicoherent supersky metric g = ^Tg“ in {na,nf,) coordinates (yellow). 


Metric ellipsoids (red) of the transformed metric Tg“, in 
max = 1500. Its projection onto the disk c = 0 (blue) generates the 


mid-time on the metrics is to rotate their metric ellip¬ 
soids (i.e. the volumes that satisfy A®A^-g^-< /Imax 
for some /Imax) around their longest axis, as may be dis¬ 
cerned from Figs. |2a] - |2el The metric ellipsoid of the semi¬ 
coherent metric g , in this case, retains the same orienta¬ 
tion of the longest axis as the coherent metrics, and still 
results in a non-constant physical metric (Fig. [^. Ap¬ 
plying T to g“, however, reorients the metric ellipsoids 


along the c axis (Fig. 2g|, corresponding to the eigen¬ 
vector of the 3x3 sky block of g“ with the smallest 
eigenvalue. Applying J then projects onto the disk 
c = 0, giving the metric g. 


Note that, contrary to Eq. (281, g is not a simple av¬ 
erage of the metrics g^, i.e. 


JT 

N 


■ N 

Es 


N N 


e=i 


e=i 


JT[re 

N 


(33) 


There are two reasons why Eq. (33) does not hold. First, 


Eqs. (27) and (28) hold only if the A® A, and hence the g^, 
are expressed in a common coordinate system; each g^, 
however, is expressed in a per-segment coordinate sys¬ 
tem (nOf,Second, the operator in general 
projects each Tg“ onto distinct hyperplanes in the coor¬ 
dinate space of n, and is therefore a nonlinear operator 
which does not commute with the averaging in Eq. (281. 


B. Validation using numerical simulations 


The ability of the semicoherent supersky metric g to 
predict the /^-statistic mismatch /I° is characterized using 
numerical simulations, the results of which are presented 
in this section. 
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1. Simulation procedure 

The simulation procedure closely follows that given in 
Section III of Paper II (which expands upon Appendix A 
of Paper I) and is briefly described below: 


(i) The simulations take a search setup as input, com¬ 

prising N segments with segment time-span T 
equally distributed over a total time-span T. The 
mean mid-time each setup is set to 

the reference time to = UTC 2015-01-01 00:00:00. 
Other input parameters are the maximum mis¬ 
matches of the coherent and semicoherent template 
banks, ^max and /imax respectively, and the fre¬ 
quency of simulated signals /. When computing 
metrics and simulating gravitational-wave signals, 
the LIGO Livingston detector [5] is assumed. 

(ii) A signal parameter vector A® = is gen¬ 

erated at reference time toi with / hxed and other 
parameters taking random values uniformly dis¬ 
tributed over the sky |n| = 1 and in spindown 
Then, the nearest semicoherent template to A'*, A, is 
determined assuming a semicoherent template bank 
constructed from an A| lattice using g with maxi¬ 
mum mismatch /imax (see Section III A of Paper II 
for details). 

(hi) The metric mismatch between A** and A assuming 
no interpolation is calculated using 


/I„i = A'^A-g. A^A. 


(34) 


The corresponding J^-statistic mismatch is com¬ 
puted, using the software package LALSuite |15] . 
by simulating a gravitational-wave signal with pa¬ 
rameters A**, computing the J^-statistic in each seg¬ 
ment at templates and A to determine the signal- 

to-noise ratios SNR^(A’*; A®) and SNR^(A’*;A), and 
calculating using Eqs. @ and ( |^ . 

(iv) In each segment the nearest coherent template to 
A, Af(A), is determined assuming a coherent tem¬ 
plate bank constructed from an lattice using g^ 
with maximum mismatch /Imax, following step ([^. 

(v) The metric mismatch between A® and A with inter¬ 
polation is calculated using 


1 w_ _ 


(35) 


The corresponding J^-statistic mismatch pP is 
computed following step (|^ by simulating a 
gravitational-wave signal with parameters A®, com¬ 
puting the J^-statistic in each segment at templates 
A® and A^(A) to determine the signal-to-noise ratios 

SNR^ (A®; A®) and SNR^ (A®; Af(A)), and calculating 


The input parameters to the simulation procedure 
are: 72 search setups given in Table |ll] two val¬ 
ues of //Hz G {100,1000}; and five combinations of 
maximum template-bank mismatches (/Imax, Mmax) € 
1(0.1,0.1), (0.3,0.3), (0.5,0.5), (0.1,0.5), (0.5,0.1)}. A to¬ 
tal of 3.1 X10® coherent J^-statistic values were computed. 


£. Simulation results: relative errors 


Following Paper 1 two mismatches ptP and /i are com¬ 
pared by plotting their relative error 


e(/r°,/r) = 


h- - h- 
0.5(/i0-k/r) 


/r®, /r > 0 . 


(36) 


If £(/r®,/i) < 0, then the metric mismatch /x over- 
estimates the true J^-statistic mismatch /x®. This im¬ 
plies that the metric is a conservative estimator of the T- 
statistic mismatch, and therefore will restrict /x® < /Xmax 
to less than the maximum metric mismatch /Xmax pre¬ 
scribed by the template bank. If e(/x®,/x) > 0, then /x 
wnder-estimates /x®; this implies that the metric is an op¬ 
timistic estimator of the J^-statistic mismatch, and there¬ 
fore may allow /x® > /Xmax- The behavior of the relative 
error can depend on how the parame ter-space offsets, e.g. 


A®A, are generated; see Figure I of Paper I If, for ex¬ 
ample, the ellipsoid A®A -g {A.^ A®) • A®A < /Xmax (where 
g°(x4.. A®) = differs in orientation from 

the ellipsoid A®A • g • A®A < /Xmax, from which A®A are 
sampled, this will tend to e(/x®;,/Xni) < 0. 

Figures [^and[^ plot the relative error between /x®; and 
/Ini, e(/IniiMni), and between /x® and /x, £(/x®,/x). Over¬ 
all, the metric g is a good predictor of the J^-statistic 
mismatch: with no interpolation, the average median 
|e(/x®i, /Xni)| ^ 0.20, and the average 25th-75th percentile 
range ^ 0.04; with interpolation, the average median 
|e(/x®,/x)| ~ 0.35, and the average 25th-75th percentile 
range ~ 0.03. The magnitudes of the relative errors with 
interpolation are generally larger than those without in¬ 
terpolation, due to the mismatch between A** and Af(A) 
being generally larger than that between A and A^(A); 
it was observed in [Paper I| that the predictions of the 
metric worsen at higher mismatches. The relative error 
magnitudes are mostly indifferent to changes in the var¬ 
ious simulation parameters plotted in Fig. 

In general, the metric g over-estimates the J^-statistic 
mismatch. The dominant contribution to this effect is the 
Taylor expansion in Fq.([^ of the J^-statistic mismatch, 
used in deriving the metric g/(^. A®) of which g^ is an 


using Fqs. ([^ and 


approximation; since g is a function of g^ by Fq. (331, 
it assumes the same approximation. Since this expan¬ 
sion is of second order, the metric mismatch /x increases 
quadratically (see e.g. Figs. 11-12 in [ID]), while the iF- 
statistic mismatch /I® is known to increase at a slower 
rate (see e.g. Fig. 10 in [20], and Fig. 7 in Paper I|. In 
addition, /x is unbounded, while /x® is by construction re¬ 
stricted to /X® < 1. This general behavior is common to 
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FIG. 3. Relative errors ef/tSi. Mni) (left column) and eijjf',11) (right column), at constant (a^ 'b) total time-span T, (c] d) 


seg¬ 


ment time-span T, and (e’^ number of segments N. All other simulation parameters are averaged over. Plotted are the 


median (solid line), the 25th-75th percentile range (error bars), and the 2.5th (lower, short-dashed line) and 97.5th (upper, 
long-dashed line) percentiles. 
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FIG. 4. Relative errors (left column), and e{'pP,'p,) (right column), at constant (a' 

template bank mismatch )imax, (c][d)| maximum coherent template bank mismatch /imax, and Je 


maximum semicoherent 
segment duty cycle 77 = 


NT/T. All other simulation parameters are averaged over. Plotted are the median (solid line), the 25th-75th percentile range 
(error bars), and the 2.5th (lower, short-dashed line) and 97.5th (upper, long-dashed line) percentiles. 
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FIG. 5. Relative errors £(/r°i, /Inj) (left column) and e{'pP,^1) (right column), at constant (a^ 'b) total time-span T, (c] d) 


seg¬ 


ment time-span T, and (e' f) number of segments N. All other simulation parameters are averaged over, subject to the 
restriction T > 1 day. Plotted are the median (solid line), the 25th-75th percentile range (error bars), and the 2.5th (lower, 
short-dashed line) and 97.5th (upper, long-dashed line) percentiles. 
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FIG. 6. Relative errors ef^niiMni) (left column), and e{pP,'p,) (right column), at constant 


template bank mismatch /Tmax, (c][d)| maximum coherent template bank mismatch /Tmax, and 


maximum semicoherent 
segment duty cycle rj = 

NT/T. All other simulation parameters are averaged over, subject to the restriction T > 1 day. Plotted are the median (solid 
line), the 25th-75th percentile range (error bars), and the 2.5th (lower, short-dashed line) and 97.5th (upper, long-dashed line) 
percentiles. 
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TABLE II. Properties of the 72 search setups used to test the semicoherent supersky metric, and of the numerical simulations 
based on those setups, as described in Section IIIB Top half: for 36 setups indexed by number of segments N (rows) and 
segment time-span T (columns), the total time-span T and segment duty cycle rj = NTjT are given. Bottom half: for 36 setups 
indexed by T (rows) and T (columns), N and rj are given. For all setups, the log-probabilities = — logj^Qp(£(/2ni,/Ini) > x) 
and = — logiQp(£(/2°, ^) > O), a; G {0,0.5}, are also given. 


all parameter-space metrics based on a second-order Tay¬ 
lor expansion; see e.g. [lOlllll) andlPaper II In addition, 


as noted above, the behavior of the relative error itself is 
biased to reporting over-estimation when the orientation 
of the metric ellipsoid of g differs from that of g°(^, A**). 

Figures and plot the same quantities as Figs. 
and but subject to the restriction T > 1 day. This 
restriction improves the predictions of g: with no inter¬ 
polation, the average median |£(/2°;,/}ni)| ^ 0.13, and 
the average 25th-75th percentile range ~ 0.03; with in¬ 
terpolation, the average median |£(/2°,/i)| ~ 0.30, and 
the average 25th-75th percentile range ~ 0.007. Many 
of the outlying relative errors seen in Figs. and [i] are 
also removed by this restriction: compare the 2.5th per- 


Figs. and vs. Figs. and and Figs. and 
vs. Figs. [4e| and [iH 

As noted in Paper I the predictions of g^ are worse at 


T = 1 day for two reasons. First, the phase metric gf{X^) 
approximates g°(^. A®) in Eq. (|^ in part by neglecting 
the amplitude modulation of the signal of a gravitational- 
wave pulsar; since these modulations arise primarily from 
the rotation of the Earth, this approximation is worst at 
T ~ 1 day. Second, the numer ical ill-co nditionedness 
of the phase metric (see m and [Paper I| increases the 


difficulty of computing g^, from which both g^ and g are 
derived, at smaller T. 

One feature of interest is the 2.5th percentile lines in 
Eigs. 5b and [if] At this line, e(/i°,/i) ^ —2 at T ^ 
10 days and N ^ 3, but reduces to ^ —1 for 

T > 40 days and N >6. The general decrease in e(/2°, fi) 
is simply due to the larger number of segments N (since 
T is approximately oc N), which mitigates the effect of a 
large over-estimation of the J^-statistic mismatch in any 
one segment. (The momentary incr ease in £(/i°, /r) along 
the 2.5th percentile line in Fig. 5b for 30 < T/days < 60 
is explained in the next section)) 


centile lines in Figs. and vs. Figs. and in ra and [5 


The jumps in the relative errors seen in Figs. ^ ^ 
most noticeably in the 2.5th and 97.5th per¬ 


centile lines - are artifacts of plotting the 72 search se¬ 
tups with respect to N. For example, note that in Fig.|3f| 
the magnitude and spread of the relative error is notice¬ 
ably less at TV = 8 than at either N = 7 oi N = 9. 
From Table |Tlj the only setup with = 8 is that with 
T = 120 days, while the setups at = 7 and N = 9 have 
T < 49 days and T < 63 days respectively. Given that 
the metric predictions improve with T, it is expected that 
the relative error at iV = 8 is reduced relative to N = 7 
OI N = 9. 
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FIG. 7. Relative errors , ^ ni) (left y-axis) at constant 
total time-span T, as in Fig. 5a Overplotted in black is the 
orientation angle P (right y-axis) of the plane of the reduced 
supersky coordinates {na,nb). 


3. Simulation results: metric under-estimation 


It is also possible for the metric g to under-estimate 
the J^-statistic mismatch. It is important to note that, 
when performing a real search using an implementa¬ 


tion of a lattice template bank, as e.g. described in Pa- 


Iper m metric under-estimation may not necessarily lead 
to loss of coverage (i.e. areas of parameter space where 
7* > Mmax)- For example, metric under-estimation at 
parameter-space boundaries may be compensated for by 
additional boundary templates. Tests of lattice template 
banks constructed using g^ in Paper II have demon¬ 
strated that near-complete coverage is attainable using 
the reduced supersky metric. 

Table [n| lists, for each of the 72 search setups, the prob¬ 
ability that the semicoherent mismatches and (i.e. 
with and without interpolation) will under-estimate the 
respective J^-statistic mismatches. The probabilities are 
given in logarithmic form by = — Iogiop(£(/J°, m) > 

x) and = - logioP(e(MnP Mni) > a;), where x G 
{0,0.5}. We observe that the probability of under¬ 
estimation for any search setup is < 10% (i.e. > 1), 

and decreases with increasing number of segments N and 
coherent time-span T, and with the addition of interpo¬ 
lation. This is expected: increasing N reduces the po¬ 
tential effect of g under-estimating the J^-statistic mis¬ 
match in a few segments, as does the addition of inter¬ 
polation which averages the coherent and semicoherent 
mismatches; and Figs. 3a and 3b demonstrate that g bet¬ 
ter approximates the J^-statistic at larger T. 


In Figures and we can see that, at the 97.5th 
percentile line, /ini noticeably under-estimates /7F for 
30 T/days < 60. This is due to a known feature 
of the reduced supersky metric; see Section V D of IPa- 


per I As the time T spanned by the metric g increases. 


the orientation of the plane of the reduced supersky co¬ 
ordinates {ria,ni,) transitions from being perpendicular 
to the equatorial z-axis to being perpendicular to the 
ecliptic Z axis; in other words, the (jia,rib) transition 
from “equatorial-like” coordinates to “ecliptic-like” co¬ 
ordinates. During this transition, the reduced supersky 
metric is more likely to under-estimate the J^-statistic 
mismatch; compare Fig. [Sa] and [5a] to Fig. 14 of jPaperTj 
To illustrate the transition. Figure [Tjoverplot s Fig. 5a] 
with the orientation angle /3 (given by Eq. (51) of jPaperT ) 
of the plane of the reduced supersky coordinates {Tia,nb), 
normalized such that /3 = 0 signifies an “equatorial-like” 
orientation and ,5 = 1 an “ecliptic-like” orientation. The 
transition from “equatorial-like” to “ecliptic-like” coordi¬ 
nates occurs between T ~ 30 days and 60 days, the same 
domain where /ini noticeably under-estimates /tF. The 
momentary inc rease in e{pP, /7) along the 2.5th percentile 
also occurs when 30 < T/days < 60, and 


5b 


line in Fig. 

therefore is also likely due to the transition. 

To determine whether metric under-estimation is cor¬ 
related with particular regions of parameter space, we 
identify pulsar sky positions, in terms of right ascension 
a and declination i5, where the metric g under-estimates 
the J^-statistic mismatch. Figures and map the prob¬ 
ability of metric underestimation as a function of (a, 5), 
for each setup given in Table [^ The maps are normal¬ 
ized such that the total probability in each ma p is either 
10“‘^ni without interpolation (Figs. 8a and 9a), or 10“'^ 
with interpolation (Figs. 


8b 


^and 9b I, where 5°; and 5° 
are given for each setup in Table^ Two effects, which 
cause metric under-estimation to be concentrated in par¬ 
ticular regions of the sky, are identified and explained in 
the following paragraphs. 

First, for short segment time-spans T < 3 days, regions 
of the sky which (anti-)correlate with the position of the 
detector used in the numerical simulations - the LIGO 
Livingston detector - are more likely to result in metric 
under-estimation; see Figure This is likely due to the 
neglect of the amplitude modulation of the gravitational- 
wave pulsar signal in deriving the phase metric g^(A’^), 
as noted above in Section [m B 2| Without interpola¬ 
tion (Fig. [8a] ), where only the semicoherent metric g is 
involved, the effect appears at T ~ 1 day and gener¬ 
ally fades for longer T and more A^; with interpolation 
(Fig. 8b), where both g and the coherent metrics g^ are 
involved, the effect is most noticeable at T ~ 3 days. It is 
probable that the effect of neglecting the amplitude mod¬ 
ulation becomes irrelevant more quickly for the semico¬ 
herent metric as N increases and more coherent metrics 
are averaged together. 


Second, at longer T, the probability of metric under¬ 
estimation becomes concentrated at either the eclip¬ 
tic equator (without interpolation), or at the equato¬ 
rial equator J = 0 (without interpolation); see Fig¬ 
ures [^ and [^ This is likely due to the transition of 
the reduced supersky coordinates from “equatorial-like” 
coordinates to “ecliptic-like” coordinates, which begins 
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FIG. 8. Intensity maps of the probabilities (a) p(e(/2ni,/2ni) and (b) p(e(^°,p) > O), as functions of right ascension a 


and declination <5, for the 36 setups listed in the top half of TableO The remaining simulation parameters, / and (/i 
are averaged over. Overplotted are the ecliptic equator (black lin^, and the positions of the LIGO Livingston detector at the 
start times of the first (black dot) and last (black cross) segments. Each map is titled by the triplet {N, T / days, r;). The 


maps are normalized such that the total probability in each map is either (a) 
for each setup in Table [H] 
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FIG. 9. Intensity maps of the probabilities (a) p(e(/2°i,/Ini) > O) and (b) > O), as functions of right ascension a and 

declination 5, for the 36 setups listed in the bottom half of Table|^ ihe remaining simulation parameters, / and (/Imax, Mmax), 
are averaged over. Overplotted is the ecliptic equator (black line). Each map is titled by the triplet {T / days, T / days, ry). 


The maps are normalized such that the total probability in each map is either (a) 
given for each setup in Table [III 
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(b) 10 ® , where 5°] and S'® are 
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at T > 30 days for the semicoherent metric (see Sec- 
and Fig. [^, and at T >10 day s for the co- 


tion 


IIIB2 


herent metric (see Section V D of Paper I|. In Figs. 8a 


and 9a (i.e. without interpolation), only the semicoherent 
metric g is in use, and the maps are generally consistent 
with a transition from “equatorial-like” to “ecliptic-like” 
coordinates on the timescale T > 30 days. In Figs. 


8b 


and 9b (i.e. with interpolation), g contributes only once 
to the total metric mismatch, whereas the coherent met¬ 
rics g^ contribute IV ^ 1 times; see the vector summation 
in Fig. [l] The maps therefore tend to reflect the behav¬ 


ior of the coherent metric, which for T <7 days remains 
“equatorial-like”. 


4- Simulation results: mismatch histograms 

Figure plots histograms of the coherent mismatches 
in each segment, Jij, Jii, and the semicoherent mis¬ 
matches with and without interpolation, fl and /i°;, 
/ini respectively. The coherent mismatches reproduce the 
results obtained in Fig. 4 of [Paper II[ The histograms 
again demonstrate the conservative over-estimation of 
the metric approximation, in that the metric mismatch 


histograms (dashed lines in Fig. 101 are peaked at higher 


mismatches than the J^-statistic mismatch histograms 
(solid lines). 

The shape of the histograms in Jli and /Ini in Figure 
conform to that expected for four-dimensional template 
banks constructed from an lattice; see Fig. 1 of jPa-j 
per II The shape of the histograms in /I, to which 


both the coherent and semicoherent template bank mis¬ 
matches contribute, is closer to a Gaussian in shape; this 
is likely due to the central limit theorem applying when 
the N coherent mismatches /I”(A;A^(A)) are averaged 
in Eq. (29), assuming N is sufficiently large. Note that, 
howe ver, t he histogram in /I for (^max, Mmax) = (0.1,0.5) 
(Fig. lOd) instead more closely resembles that of /Ini, as 


might be expected since flmax 'C /Imax, he. the coher¬ 
ent template bank mismatches fig contribute little to the 
total mismatch /I. 


IV. COMPARISONS TO PREVIOUS WORK 


Section IV A compares the semicoherent metric g, pre¬ 


sented in the previous section, to the global correlation 
metric of [28] . Section IV B compares predictions for the 


number of semicoherent templates given by g to predic¬ 
tions given by UHl HH] • 


A. Comparison to the global correlation metric 

The numerical simulations described in Section IIIIBI 
were repeated with two changes: the coherent global cor¬ 
relation metrics gf'^ given by Eqs. (50) of [55] were substi¬ 


tuted for the coherent supersky metrics g^, and the semi¬ 
coherent global correlation metric g®'^ given by Eqs. (59) 
of |55| was substituted for the semicoherent supersky 
metric g. In other words, the simulations were repeated 
assuming that the global correlation metrics are used to 
construct the coherent and semicoherent template banks 
in steps (3 and (|^ of the procedure described in Sec¬ 
tion [111^ 

Figure ^ plots the relative error between /I°j and 


and between JjP and /I®“, e{fjr, /I®“). Figs. 


11a 


and jllbj are directly comparable to Figs. and [3b for 
the supersky metrics respectively. The global correlation 
metric over-estimates the J^-statistic mismatch, with the 
median |£(/Ini,/Ini)l ~ |£(m°,^ 0-4, provided that 
T < 20 days; for T > 20 days, the global correlation met¬ 
ric severely under-estimates the J^-statistic mismatch, 
with the median £(/lF,/I®;) > 1 and £(/I°,/ls“) > 0.5. 
Figs, jllcljllfl restrict the search setups listed in Table [III 
to T e {1,3} days and r/ = 1 (i.e. contiguous segments). 
The behavior of the median relat ive error is similar ^n- 
der these restrictions: in Fig. lld[ £(/I°,/I®“) > 0 for 
20 days and £(/I°, /!«“) > 0.5 for T > 120 days; in Fig. 

> 0 for T > 14 days and £(/I°,/I®“) > 0.5 for 
T > 30 days. 

The relationship between the semicoherent and coher¬ 
ent metrics is often quantified using the refinement, 


Ilf 


7 

This quantifies the ratio of the number of templates in 
the semicoherent template bank Af oc -^/detg to the me¬ 
dian 1451 number of templates in the coherent template 


plots the refinement for both supersky and global corre¬ 
lation metrics, the former given by Eq. ( [37[ ), the latter 
given by Eqs. (50), (^, and (73) of |28|. The search 
setups listed in Table |ll| are restricted to T = 1 day and 
rj = 1. For T > 20 days, the rate of increase of the su¬ 
persky metric refinement increases markedly, while that 
of the global correlation metric refinement remains con¬ 
stant. The lack of additional refinement by the global 
correlation metric for T > 20 days is consistent with the 
deterioration in its ability to predict the J^-statistic mis¬ 
match above the same threshold, as observed in Fig. El 
The semicoherent global correlation metric is derived 
under a number of assumptions, which were studied in 
unpublished work by |46| . One assumption is that the 
sky position parameters do not contribute to the refine¬ 
ment, i.e. that the number of distinct sky positions that 
must be searched remains the same in both the semico¬ 
herent and coherent template banks. We conclude from 
Figs. [^ and [^ however, that this assumption remains 
valid only for T < 20 days. That the sky position param¬ 
eters do contribute to the refinement has also been found 
in work on the Hough method, a similar semicoherent 
search method for gravitational-wave pulsars [IZliTj. 


banks A/} oc y^det g^, assuming ^max = Atmax- Figure 
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V^det) 


median^} det g^} 


(37) 
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^4, fit 



/^ni’ /^ni 


(a) 



fJ-i 



/^ni’ /^ni 

(b) 



ttO t; 
f^i 1 



/^ni’ /^ni 

(c) 




^ii /^ni’ /^ni 
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—0 t; 
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FIG. 10. Histograms in mismatches in: the coherent mismatches in each segment Jli (left column); the semicoherent 
mismatches with no interpolation /r®;, /Ini (middle column); and the semicoherent mismatches with interpolation /t°, Jl (right 
column). Hist ograms a re a t fixed maximum template-bank mismatches (/Imax, Mmax) equal to [R] (0.1,0.1), [^(0.3, 0.3) 


(0.5,0.5), (d) (0.1, 0.5), (e) (0.5,0.1). All other simulation parameters are averaged over. 
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T / days T / days 

(e) (f) 


Relative errors e(/rni!A‘S) (1®^^ column) and (right column), at constant total time-span T. All oth er 

~ l(a][b) 


FIG. 11. 

simulation parameters are averaged over, subject to the following restrictions on the search setups listed in Table II 


restrictions; (c’ d) T = 1 day, ?? = 1; and (e' f) T = 3 days, 77 = 1. Plotted are the median (solid line), the 25th-75th percentile 


range (error bars), and the 2.5th (lower, short-dashed line) and 97.5th (upper, long-dashed line) percentiles. 
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/max5 'min5 


given by g^. Here, using the same values of /, 
and //max = //max, the p rediction for Af given by g, which 
follows from Eq. (29) of Paper III is compared to the pre¬ 
diction of [55], given by Eqs. {2.15)-(2.19) and (2.22)- 
(2.27) of that paper, and to the prediction of the global 
correlation metric, given by Eqs. (56), (59), and (65) 
of |5B|. The search setups listed in Table jnj are restricted 
to r = 1 day and rj = 1. 

In |26j . the number of semicoherent sky templates is 
computed assuming that only the rotational motion of 
the Earth is important; a similar assumption is made 
in deriving the global correlation metric |28j . Paper I 


T / days 


demonstrated, however, that both the rotational and or¬ 
bital motion of the Earth are important to the coherent 
metric, and the same is true of the semicoherent metric. 
The neglect of the orbital motion by both [5^ and |5S] is 


FIG. 12. Refinement of the supersky metric given by Eq. ( |37[ ) 
(solid line), and of the global correlation metric given by 1281 
(dashed line), as a function of total time-span T. The search 
setups listed in Table [n| are restricted to T = 1 day and 77 = 1. 



T / days 

FIG. 13. Number of semicoherent templates predicted by the 
supersky metric (solid line), by the global correlation metric 
of |28| (dashed line), and by [55] (dash-dotted line), as a func¬ 
tion of total time-span T. The search setups listed in Table [n| 
are restricted to T = 1 day and 77 = 1 . 


B. Number of semicoherent templates 


consistent with their predictions for Af not following the 
increase predicted by g for T > 20 days, once the orbital 
motion begins to dominate, as seen in Fig. |13[ 

In Section IIIIB 21 we saw that the semicoherent met¬ 
ric g generally over-estimates the J^-statistic mismatch 
pP. The metric will therefore generate a template bank 
with more templates than are strictly required to sat¬ 
isfy the condition pP < ^max- This effect cannot, how¬ 
ever, account for the increase in the number of templates 
seen in Figure |13[ for two reasons. First, the number 
of templates scales with the maximum metric mismatch 
as Af (X ^max^, where n = 4 for a four-dimensional 
parameter space. Let us assume very conservatively 
that g always over-estimates JA’ by a factor of 2, i.e. 

= -0.667 (cf. Figs. I and|^. We could there¬ 
fore increase //max by a factor of 2 to achieve the desired 


maximum /z*^, which would reduce the number of tem¬ 
plates by a factor (//max/2)“^/^///max^ = 4. This factor 
is still much smaller than the 10~^ increase in the num¬ 
ber of templates predicted by g, seen in Fig. ^ this 
increase cannot, therefore, be accounted for by metric 
over-estimation. Second, as discussed in Section jlll B 2 


the dominant contribution to the over-estimation of //' 
by g is the Taylor expansion in Eq. ([^ used in deriving 
the phase metric. Both the metrics of [26] and of [28], 
against which g is compared in Fig. [^ are also derived 
from a phase metric via Taylor expansion. The met¬ 
rics of [26l I 2 H] will therefore, in addition to their other 
assumptions, also tend to over-estimate the J^-statistic 
mismatch. 


Figure [T^ plots the number of semicoherent templates 
Af required to search the whole sky, frequencies up to 
/max, and spindowns |/| < //Rnin (motivated by a min¬ 
imum “spindown age” of the gravitational-wave pulsar), 
with a prescribed maximum mismatch //max- This pa¬ 
rameter space is used in [15] and [55] for predictions of 
the number of coherent and semicoherent templates re¬ 
spectively. Setting /max = 200 Hz, Tmin = 10^ yr, and 
/tmax = Mmax = 0.3, the prediction for t he number of co¬ 
herent templates in m was compared in||Paper H[to that 


V. SUMMARY 

This paper presents a useful approximation to the 
semicoherent parameter-space metric: the semicoherent 
supersky metric g. Together with the coherent super¬ 
sky metric introduced in [Paper l[ it constitutes a com¬ 
plete description of the parameter-space metrics for all¬ 
sky semicoherent searches for isolated gravitational-wave 
pulsars. The number of semicoherent templates required 
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to maintain a prescribed maximum mismatch is found 
to be several orders of magnitude larger than previous 
predictions. 


The consequences of the required increase in the num¬ 
ber of semicoherent templates for the sensitivity of 
gravitational-wave pulsar searches is left for future work. 
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